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^ - Abstract. The properties of some locally rotationally symmetric (LRS) perfect fluid 

^ , space-times are examined in order to demonstrate the usage of the description of 

I geometries in terms of the Riemann tensor and a finite number of its covariant 

ly^ ' derivatives for finding solutions to Einstein's field equations. A new method is 

introduced, which makes it possible to choose the coordinates at any stage of the 
calculations. Three classes are examined, one with fluid rotation, one with spatial twist 
I in the preferred direction and the space-time homogeneous models. It is also shown 

that there are no LRS space-times with dependence on one null coordinate. Using an 
. extension of the method, we find the full metric in terms of curvature quantities for 

the first two classes. 
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PACS numbers: 04.20.-q, 04.20.Jb, 95.30. Sf, 98.80.Hw 
1. Introduction 

Searching for exact and approximate solutions to Einstein's field equations can be a 
cumbersome task, and often ad hoc mathematical assumptions are used to solve the 
obtained equations. This leads to the problem whether or not the solution is physically 
interesting, which need not be obvious. It is therefore convenient if one can start from 
physical restrictions, e.g., on a fluid say, and from that derive a solution to Einstein's 
field equations. 

In general relativity the Riemann tensor and its covariant derivatives is essentially 
what one can measure, i.e., those are the physical observables. From the equivalence 
problem |^ we know that given a metric g = rjijuj^ ® in a moving frame, the set 
RP+^ = {Rijki, Rijki-mi...mp+i}, consisting of the Riemann tensor and a finite number 
of its covariant derivatives, gives a complete local description of the manifold (p -|- 1 is 
defined as the lowest number for which the derivatives of that order are functionally 
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dependent of those of lower order). This can be used for classifying metrics, roughly 
by comparing the different sets for different metrics 0, This naturally leads to 

the question of whether or not it is possible to construct the local geometry (i.e. if we 
can find some basis and connection 1-forms satisfying Cartan's equations) of a manifold 
starting from some set of elements which we want to be the Riemann tensor and its 
covariant derivatives up to some order 0, ||, If. This task was pursued in Refs. 
and it was shown that one can always do this if some integrability conditions (parts 
of the commutators, Ricci identities, cychc and Bianchi identities) are satisfied by the 
elements of the set. 

The method has not yet been used in an extensive way for finding solutions to 
Einstein's field equations. Therefore, in this series of papers we will study how the 
method applies on locally rotationally symmetric (LRS) perfect fluids. This is quite 
a natural choice, since it has a well defined symmetry property and contains several 
physically interesting classes, e.g., inhomogeneous models and spherically symmetric 
models of astrophysical interest (to be studied in Ref. 0). Also, it has been studied by 
other means, thus making a comparison between different methods possible. 

In Sec. H we review the equivalence problem and its inverse procedure, i.e., 
construction of solutions of Einstein's equations. 

In Sec. 0we set up the conditions for a perfect fluid space-time to have LRS, and we 
generate the integrability conditions (IC) for these space-times. Here we introduce a new 
technique for keeping the choice of coordinate (s) undetermined through the calculations. 
We classify the space-times into six classes which are either disjoint, or need some special 
treatment. In Sees. |^-|^ we investigate the IC for four of the classes. 

In Sees. [7|-PD| we generate and solve the isometry algebra over F{M) for two classes, 
and find the form of the metric for these. 

We will use the following conventions: 

• The metric has signature (+,—,—,—). 

• Uppercase latin indices refers to the frame bundle basis: 

A, B, ... = 1, 2, k, 

P, Q,... = k + l,k + 2, n{n + l)/2, 

and I, J, ... are general indices. 

• Lowercase latin indices denotes components with respect to the moving frame on 
M: 

••• = 1,2, ...,n. 

Boldface lowercase latin indices refers to the functionally independent objects in 

I The method has been given the name Karlhede classification scheme (and sometimes the Cartan- 
Karlhede classification scheme). 



3 

ci, ... ^ X , 2 , . . . , k . 

• Uppercase greek indices is used for coordinates on the fibers of tlie frame bundle, 
i.e., the parameters of the orthogonal group: 

T,$,... = l,2,...,n(n-l)/2. 

• Lowercase greek indices are used for the coordinates on M: 

fi, z/, ... = 1, 2, ...n. 

• f\i denotes the directional derivative ej(/), where / is some function and ej a 
tangent basis vector. 

2. Method 

Here we first give a short review of the equivalence problem and its solution. Suppose we 
are given two tt,- dimensional Riemannian manifolds (M, g) and (M, g) . The two metrics 
are equivalent on some subsets U C M, U C M ii g|[/ = g|^. In terms of coordinates, 
the metrics become g = g^ydx^ ® dx*^ and g = g^ydx^ ® dx" , and two geometries are 
equivalent if there exists coordinate relations x^ = x'^{x'^), /x, z/ = 1, 2, n, such that 

Ox'^ Ox^ 

Since it is difficult to find these coordinate transformations in general, this approach is 
not effective for investigating equivalence. 

If one uses a moving frame approach, i.e., one introduces a tangent basis {ej} and 
its dual {w*}, i = 1,2, ...,n, (and analogously for M), so that 

g = riijuj' ®u)^ , g = r]ijCj' ® G)^ , (2.2) 

where {rjij) is some constant, symmetric and nonsingular matrix, the equivalence 
problem can be put in a new guise. From (|2.2|) we see that a>* = u)"^ implies g = g. 



Unfortunately the converse is not true. If G is the orthogonal group on M, we can 
use some representation of G to make transformations of the basis while keeping the 
tetrad components of the metric fixedf, i.e., 3 some A G G such that r]ki = ly^h.h'Uij (we 
will not bother in distinguishing between G and representations of G). This freedom 
is captured by letting the basis 1-forms be defined over the frame bundle F{M), i.e., 
they depend on both the coordinates {x^} over the manifold M and the parameters 
{^^}, T = 1,2, ...,n{n — l)/2, of G. Cartan's equations of structure, which defines the 



f Actually, there still exist the possibility of discrete transformations (spatial reflection etc.) but since 
these are finite in number they are rather trivial (at least in our case). We therefore suppress this 
freedom. 
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connection 1-forms uj^j and the curvature 2-forms R*j =^^i?*jfc;a;'^ Ao;', on a Riemannian 
manifold reads 

da;* = a;'^ A uj\, (2.3) 
du;'j = -u:\ A UJ^ j + YCj . (2.4) 

with the condition rjikOJ^j + rjjkU^'^i = 0. Since our 1-forms are defined over the frame 
bundle, the exterior derivative is defined as d'==d2;+d^, which implies that the connection 
1-forms consists of two parts: 

a;* =^'W, , (2.5) 

where ^^^cx^^j'^T'^jku:'^ and ^^-'o'^- =^ajTd.^^, for some Fj^ and a^x determined by Cartan's 
equations. Here T^jk are the Ricci rotation coefficients. Therefore, when defined over 
F{M), the connection 1-forms are linearly independent of the basis 1-forms. They make 
up the missing 1-forms for the basis {a;^}, J = 1, 2, n(?2 + 1)/2, of the cotangent space 
T4F(M)], i.e., {u;^} = {u;\uj)}. When dx'^ = 0, i.e., along a fibre in F(M), the ^^^'^uj) 
reduce to the generators t) of G and they fulfill d^r j = — A r^^. 

Using this, one can state the equivalence problem in the following way |1], Q]: Two 
geometries are equivalent if and only if 3 relations x'^ = x'^(x'^), = ^^(x'',^*) such 
that a>*(x'',f^) = uj^{x>',^'^). 

With the collective notation introduced above for the basis of T*[F(M)], we can 
write Eqs. (2.3) and (2.4) as 

da;^ = lcW^Aa;^ , (2.6) 

where C^jk essentially represents the Riemann tensor. In Refs. [|T], |^ it was proved that 
two metrics g and g are equivalent if and only if 

^ JK — ^ JK , 

ril /-(/ 

JK\Ni — .JK\Ni 1 



C^jK\Ni-Np+i - C\jK\N^-Np+^ , 

are compatible as coordinate relations over F{M). Here p + 1 is the lowest derivative 
order for which C^jK\Ni-Np+i is functionally dependent on the derivatives of order < p+1. 
This implies that the set 

^ — \^ JKi JK\Ni, • • • , O JK\Ni-Np+if 

gives a complete local description of the manifold. In terms of the Riemann tensor, this 

set is given by RP+^ = {Rijkl, Rijkl;mi, ■ ■, Rijkl;mi-mp+i} 0- 
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Since we can achieve a local description of M in terms of RP^^, it can be used for 
constructing solutions to Einstein's equations in an invariant way, because we know 
that two solutions of Einstein's equations are equivalent if a comparison between their 
^p+i gives consistent coordinate relations. Obviously, any set of RF^^ elements can not 
generate 1-forms such that Cartan's equations are fulfilled. Some conditions have to 
be imposed on the set. Let {/^}, a = l,2,...,k < n{n + l)/2, be a maximal set of 
functionally independent objects chosen from W [this means that the dimension of our 
isometry group is n{n + 1)12 — k]. As our basis we have chosen {cj^} = {cj*, ct» j}, where 
{a;*} is some moving frame. In Refs. [1], it was shown that the integrability conditions 
for a set are 



d'/^ = d (/V^ j = , (2.7) 
d^o;^ = d( \c''jK^^ ^uA=Q , (2.8) 



.2 

where P = + 1, + 2, ...,n{n + l)/2. When these are fulfilled, we can solve for {u;^}, 
A = 1,2, ■ ■ ■ , k, through Cartan's equations. In Ref. a procedure for finding the rest 
of the 1-forms, based on Ref. 0, was described. If the manifold is without symmetries, 
i.e., k = n{n + l)/2, we see that we do not need any of Eqs. (2.8). If we impose some 
symmetry, we reduce the number of functionally independent elements in RP^^, and 
Eqs. (2.8) have to be added. 

In practice it is often easier to work in a fixed frame than to let the components in 
i?^"*"^ depend explicitly on the parameters of the orthogonal group [i.e., we choose 
a cross-section of F(M)]. Suppose that R^'^^ only depend on , a = 1,2, ...,/, 
and rotations in the a6-planes, {°'b} = l,--.,m, where I = n — dim(orbits) and 
m = n{n — l)/2 — dim(isotropy group). Equation (2.7) then correspond to [^ |^ 

d^x" = d{x\u:') = ^ [e,, ei](x°) = -{T^i - Fifc)e,(x") , (2.9) 

dr^ft = d (uj"'b — r^bi'*^') R^ij = 2 e[j(r"|b|j]) + r''m[iX"^\b\j] + ^%k^\ij] , (2.10) 

where | means exclusion from the antisymmetrisation, and the rest of the equations that 
have to be satisfied are (2.8) which can be written 

R) Auj^ = ^ R\ijk] = , (2.11) 

dRP, + RPfe A - uj\ A R^, = ^ R%,,u] = , (2.12) 

where t = 1 + 1,1 + 2, ...,n and {^q} = m + 1, m + 2, ...,n{n — l)/2. For general relativity, 
the usual 1+3 orthonormal frame approach (i.e. splitting of space-time with respect to 
a timelike congruence) often makes use of the field equations together with the Jacobi 
identity for the tangent basis. Here they are replaced by Eqs. (2.10) and (2.11) (we 
assume that all other symmetry properties of the Riemann tensor are fulfilled). 

In Ref. [^ it was also shown that instead of using the full set R^^^, one can use 
the reduced set S = {R^jki,T%k,x^\k,riij} to describe the geometry, since PP"*"^ can be 
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constructed from 5*. Note the difference in that the set R^^^ is covariantly defined, 
while the set S is not (at least in the generic case). 

From the results listed above we see that if the manifold have translational 
symmetries only, we do not need any of the Bianchi identities. When we add some 
kind of isotropy we need some of the Bianchi identities but we can reduce the number 
of Ricci identities needed. 

Thus, to construct a geometry given a set R^^^ or S, we adopt the following scheme: 

1. Determine the set R^^^ or S, with the help of the IC. 

2. Determine the uj^. Since we can write dJ^ = I^\iuj^ , we can invert this relation for 
k of the 1-forms, i.e., for {u;^}. Since the inverse of are part of R^~^^, these 
functions are determined from the IC. 

3. Derive the isometry group from the projected Cartan's equations. 

4. Finally, determine uj^ as follows: 

(a) Using the isometry group, solve for the 1-forms in terms of some coordinate 
basis. Note that we can always pose this as a boundary value problem for a set 
of coupled ordinary differential equations (see the constructive proof of Lie's 



Third Theorem in Flanders [|T0|)- In most cases, there already exists canonical 



choices tabulated, e.g., in the 3-dimensional case the Bianchi classification [11 



and in the 4-dimensional case the work of MacCallum |]T2[. Now, use one of the 
procedures above to make a metric ansatz g by extending the 1-forms found 
to the entire manifold, 
(b) From this, calculate the set R^~^^ and compare with the original set, giving 
consistency equations for the coefficients in the metric ansatz. 



3. Integrability conditions 

LRS perfect fluids have been extensively studied, because they are simple in their 
symmetry and contains a lot of physically interesting examples (see Refs. [[ 

21). 



16, 17, IS, ic,pq, 21 
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The energy-momentum tensor is given by 

T = (yu + p)u ®u-pg , 

where u is the 4-velocity of the fluid, /i the energy density, p the pressure and g the 
metric. We choose a comoving Lorentz-tetrad, i.e. we select the tangent basis {ej} and 
its dual {uJ^}, i = 0, 3, such that 



.2\2 



3\2 



U 



(3.1) 
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We have the following definition of LRS [0, 0, ^ : 

Definition A space-time is said to be LRS in a neighborhood U{p) of a point p & M if 

at every point q G U{p) there is a subgroup of the proper Lorentz group that leaves in- 
variant the Riemann tensor and its covariant derivatives, i.e., there exists a continuous 
isotropy group at each point in the neighborhood U (p) . 

We rotate the axes so that our symmetry lies in the 23-plane, which makes our 
tetrad fixed up to rotations in that plane. 

Assumptions: The essential coordinates on the space-time manifold M are x*^ and 
x^, where x^ is time- like and x^ is space- like, f We assume that their exterior derivative 
can be written^ 

dx° = Xa;° + (3.2) 

dx^ = xuj° + yu}\ (3.3) 

where X = x°|o, Y = a;°|i, x = x^\o and y = x^\i. If the coordinates are independent we 
must have Xy — xY ^ 0. Since is time-like, i.e., X"^ — > 0, we must have X 7^ 0. 
An analogous argument gives that y ^ 0. From this we also see that if X = or ?/ = 0, 
then y = or X = respectively. 

In previous papers 0, ||, ||, the coordinate have been specified in advance (for ex- 
ample, the density). A drawback with this procedure has of course been that one has 
then restricted once attention to models with some special behavior (in the previous 
example of coordinate, we must not deal with an incompressible fluid). This drawback 
is eliminated by using the general structure introduced above. Also, since there often 
exists several choices of coordinates in a given situation, this will enable one to make 
the 'canonical' choice, i.e., the choice that gives the simplest equations. 
The reduced set S contains the Riemann tensor 

i?oioi = E - + 3p)/6 , i?0202 = -R0303 = -E/2 - (/U + 3p)/6 , 

-R1212 = -R1313 = E/2 — , -R2323 = —E — /i/3 , 

-R0123 = 2-Ro213 = — 2i?0312 = H , 

where E = En and H = Hn are components of 

E^^Ck^iu^'v' , (3.4) 
H,,=e,umC''^jnu''u''/2 , (3.5) 

t These are elements of , or combinations thereof. 

X Since the a;"^ are curvature quantities from the set RP^^ , components in the 62 and 63 directions 
would break the LRS. 
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where Eij and Hy are the electric and magnetic parts of the Weyl tensor respectively 
(here eikim is the totally antisymmetric tensor defined by €0123 = !)• 

To be able to construct a geometry we also need the rotation-coefficients. Some of 
these are expressible in the kinematic quantities acceleration a^, expansion O, vorticity 
Uij and shear aij through 

—Toij = VjUi = uJij + aij — hijQ/3 + aiUj. (3.6) 

The kinematic quantities are defined as |^ 



©def I— 7 ?■ def 



and hij'^UiUj — rjij is the projection tensor. The choice of a comoving frame results 
in all kinematic quantities with a index being zero, = uJiQ = cTjo, because of their 
orthogonality to the 4-velocity. Because of the LRS there is only one independent 
component of the vorticity, UJ23 = —1^32 = ^, and one component of the acceleration, 
ai = a. The shear will, because it is traceless, look like (0"^) = diag(0, — 2(T, a, (j).The 
rotation coefficients are then 

Toio = -a , Foil = 2a + 0/3 =^ a, 

r022 = r033 = + 0/3 = ro23 = ~r032 = 5 

-p -p def -p p def \ 

J- 122 — J- 133 — ~ ; J- 123 — 132 — — A , 

where the restriction on the ri2j and Ti-^i follows from the LRS requirement. We do not 
need r23j since we have rotational symmetry in the 23-plane. The rotation coefficient k 
corresponds to the spatial divergence of the vector field ei, while A is the spatial rotation 
of the same vector field relative to a triad which is Fermi-transported along u. Thus all 
rotation coefficients are covariantly defined because of the LRS. 

We can observe that we have a 3-dimensional spatial isotropy group if and only 
if E = H = a = uj = a = y = x = Y = 0. This also means that we do not need the 
ri2i and Fisj, and some of the Ricci identities are redundant. On the other hand, we 
need further components of the Bianchi identities (these will turn out to give no more 
information). Spherically symmetric models must have LRS and an acting only at 
one point. Therefore, at that point (the 'center' of the model) all quantities listed above 
must be zero. They then evolve smoothly from the center, and the evolution is given 
by the IC. 

With all this we are ready to let the IC take explicit form, and we give them in the 
Appendix. It is interesting to note that to determine the magnetic part of the Weyl 
tensor, H, we do not need any evolution or divergence equations. It is algebraically 
related to the other quantities in [see Eq. (A6)]. This is because we do not need 
the components of the Bianchi identities where H appear differentiated (according to 
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the theorem given in Refs. f^, 0). In general, space-times without any isotropy have 
their Weyl tensors algebraically determined in terms of the rotation coefficients, their 
derivatives and the coordinate gradients through the Ricci identities. This follows as a 
corollary of the theorem presented in Refs. 0] . 

Although they are not needed (since they are contained in the IC), things will be 
greatly simplified if one uses the twice contracted Bianchi-identities: 

Ai|o = -(yu + p)0, P|i = (/^ + p)a- (3.7) 

These can be obtained by applying eo and ei on parts of the Ricci identities, and using 
the rest of the IC. 

We can compare the curvature description of geometries used here with the 1+3 
threading formalism presented in the paper by van Elst and Ellis ||2^ , since their article 
resembles the ffist part of this one (for a general discussion of the 1+3 orthonormal frame 
approach, see Ref. [24| ). Using their notation, they choose a vector as their preferred 



(normalized) space-like vector field, and defines the rotation of it in the usual manner. 
Since the rotation has to be proportional to e*, they have e*-''^' (VjCfc) ui = —ke^ , for 
some function k. If we translate this to our notation, i.e., choosing a comoving Lorentz 
tetrad, and using e* = 5\ (using their signature convention) we get that k = 2T132 = 2A. 
Further, the spatial divergence of e\ avEE'= ^^iVjC-^ (index by the author), will satisfy 
OvEE = 2k. The procedure of obtaining the relevant equations differs somewhat between 
the 1+3 formalism and the curvature description. In the former, one generate the 
Ricci and Bianchi identities together with the Jacobi identities. Independence of the 
equations is then checked. One obtains a split in the above equations into the evolution 
equations and the constraint equations (which contains only spatial derivatives). From 
this, one checks the consistency by taking the covariant time derivative of the constraint 
equations and demand that this derivative vanish. In the Bianchi identities, there are 
equations containing the derivatives of H. Now, we can from the curvature description 
conclude that these equations are redundant, since (because of the LRS) these parts of 
the Bianchi identities are never needed. Although this is is no big advantage in this case 
(because in the three cases below the Weyl tensor is either algebraically determined 
or has zero magnetic part) it points at the fact that the reduction of the number of 
equations is automatic. 

A classification of the different LRS models can be done using the IC. One essentially 
uses Eqs. (A2) and (A4), together with the causal properties of the coordinates (i.e., 
that is time-like and is space-like). 

First of all, if the rotation is nonzero we can use Eqs. (A2) and (A4) to write 

X = YX/uj and x = yX/uj (3.8) 

respectively. Inserting this into the causal restrictions on the coordinate gradients, we 
see that the equations sphts into two cases (i) x = y = or (ii) X = Y = (there 
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is also the null coordinate case — = — = 0, which we treat later). Case 
(i), using Eq. ( |3.8| ) to combine Eqs. (A9) and (AlO), gives /? = —k\/uj. Inserting this 
into Eqs. (A15) and (A16) gives (/i + p)k\ = 0. Discarding the case /i + p = for 
the moment, we check k = 0, and see that A = 0, so A = always holds [if k = 0, 
Eq. (All) gives Xu = 0]. But this gives X = 0, which in turn implies Y = 0, i.e. no 
dependence on any coordinates at all. These space-time homogeneous cases is treated 
in Section ^. Performing an analogous manipulation for (ii) gives x = 0. But this does 
not imply y = 0. Thus we here have a nontrivial case. From Eq. ( |3.8D we get A = 0. 
Inserting this into Eqs. (A9) and (A14) results in /5 = a = 0, i.e., the fluid is shear- and 
expansion- free. 

Second, we can have A 7^ 0. From (A2) and (A4) we have 

Y = Xuj/X and y = xuj/X, (3.9) 

respectively. If we insert this into the causal restrictions on the coordinate gradients, 
we obtain (as in the first class) two separate cases: (i) x = ?/ = or (ii) X = Y = 0. For 
(i), using the relation between X and Y, we can combine Eqs. (A13) and (A14) and get 
K = —j3uj/X. Making use of this, Eqs. (A15) and (A16) gives the result (yU + p)j3uj = 0. 
As before, we assume + p ^ for the moment. We see that this implies that k = 0, 
and multiplying Eq. (All) with u gives u = 0, i.e. Y = 0. From Eq. (A9) we obtain 
a = 0, which also can be seen from Eqs. ( p.7| ). We can proceed in the same manner 
for (ii), and obtain y = 0. But this implies a; = (see Assumptions), i.e., space-time 
homogeneous models (see Section ||). 

Third, if a; = = A we can have dependence on both a time-like and a space-like 
coordinate in general. This case have zero magnetic part of the Weyl tensor, which can 
be seen from Eq. (A6). 

When fi + p = 0, the IC implies that fi and p are constants. Therefore this is 
equivalent to vacuum with a cosmological constant, and we will from now on assume 
that /i + p 7^ 0. 

Thus, the IC gives (when fi + p ^ 0): 

1. 7^ (LRS class I) ^X = F = x = A = cr = = O, no time-like dependence. 

2. A 7^ (LRS class III) =>x = y = Y = K, = a = uj = 0,'n.o space-like dependence. 

3. a; = = A (LRS class II) ^ if = 0, generally both space-like and time-like 
dependence. 



The classification in the parenthesis is the one given in Refs. [[13], |T4|, |2^. Three cases 
which lies somewhat outside the above classification are 

4. Space-time homogeneous cases. 



5. Dependence on one null coordinate. 
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6. yU + p = 0. 

The fourth class contains e.g. the Godel universe and the fifth class will be shown to be 
homogeneous. 

In this paper we study the classes 1, 2, 4 and 5. The third class contains e.g. 
the spherical symmetric case, which means that there is an abundance of models in it. 
Therefore we treat it separately in a forthcoming paper 

4. Nonzero vorticity: a; 7^ 

In this case, we see that many of the kinematic quantities are automatically zero. 
This of course means that our system of equations is heavily reduced. It in fact 
becomes two parts, one consisting of differential equations, and one consisting of 
algebraic relations. Inserting the algebraic relations into the differential ones, we get 



the differential equations 

yuj' = (2k, — a)uj, (4-1) 

yK = (fi + p) /2 ~ an — Lu'^ + k^, (4.2) 

ya' = - + 3p) /2 + a^ + 2aK + 2uj'^, (4.3) 

yp' = a{fx + p) , (4.4) 

where ' denotes differentiation with respect to x^, and the algebraic relations (which can 
be seen as defining relations) 

E = - (/i + 3p)/3 + 2aK + 2w^ (4.5) 

H = 2{a-K)u. (4.6) 



Equations (4.1)-(4.4) contain a set of six functions: {u, k, a, y,p, /i}. One of them (or 
a combination, if preferable) can be chosen as coordinate, which makes y algebraically 
determined in terms of the other functions. If we assume an equation of state the number 
of functions reduces to 3 (counting only those contained in the differential equations), 
i.e., the same as the number of equations. Thus, in the generic case, we need only solve 
three differential equations, and the rest of the functions are determined algebraically. 

If there are some algebraic constraints (e.g. E = 0) we can insert this in the IC. If 
no equations turns out to be linearly dependent, and we assume an equation of state, 
we can use this to generate a second order differential equation for the equation of state, 
as will be seen below. All other quantities are then algebraically determined in terms 
of the pressure and density (and possible constants), when choosing the pressure as x^. 

We can observe that for conformally flat space-times, the constraints on the 
kinematic quantities gives (when inserted into the IC) as equation of state fi + p = 0. 
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4.1. Rigid rotation 

If we have rigid rotation, i.e. u' = 0, Eq. (4.1) gives a = 2k,. This can be used to 
combine Eqs. (4.2) and (4.3) into an algebraic equation for k^: 

k'^ = {3fi + 5p)/20-2uj^/5 , (4.7) 

which, when differentiated (assuming an equation of state), yields (i) /t = 0, or (ii) 
2/i — 12^"^ — 3{fi + p)fi,p = 0. Case (i) will be treated for general vorticity in Sec. |4^ . 
For case (ii), the (inverse) generic solution becomes 

p{fi) = 12uj^ + A\6iu^ - fxf^ , (4.8) 

where A is a constant of integration. There is also the solution /i = Qu"^. In general, we 
may choose = p, which gives y = 2k{ii + p). 



4-2. Space-times with constant pressure 

From (4.4) we see that a = 0. From Eq. (4.3) we then get a;^ = {fi + 3j9)/4. 
Differentiating this expression for u and using the IC, we get 

yfi' = 4K(/i + 3p) . (4.9) 

If the energy density is constant, we get a homogeneous model (see Sec. |^). 
If the energy density is non-constant, we can choose this as our spatial coordinate. 
From (|4^ ) we see that y = 4/t(/i + 3p). We then obtain 

K^fx) = (/i + 3p)/2 + 2p + A{fx + 3j9)^/2 ^ 

E(/i) = {i2 + 3p)/Q , H\fi) = (/i + 3p)V2 + A{fi + 3pf'^ + 2p{^i + 3p) . 



Comparing with van Elst and Ellis |2^ for dust space-times, we see that it possible to 
integrate the system with the above choice of coordinate, thus eliminating the procedure 
of solving the second order equation (67) in [p^] for /z. 



4-3. Vanishing electric part of the Weyl tensor 

= is equivalent to the constraint cu^ — (/i + 3p)/Q + an = 0. Differentiating this (and 
using the IC) gives a constraint on the equation of state /x = /i(p): 

a(yU + p)/i,p = 3(k — a){fi + 3p — Gan) . (4-10) 

This equation can be used to replace (4.1), so that uj is determined hj E = 0. From 
this we see that a constant energy density gives a = k, thus making the space-time 
conformally flat. But this implies u; = 0, contrary to our assumptions. 

The general system one obtains is fairly complicated, since it contains 9*'' degree 
polynomials in a and k, with coefficients consisting of p, /i and its derivatives up to 
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third order with respect to p. It therefore seems unhkely that it is possible to find 
the general equation of state. A possible solution to this problem is to introduce some 
further constraint, that will simplify the calculations. We might insert a linear barotropic 
equation of state p = (7 — The procedure above then gives a polynomial in 7 — 1, 
in which all coefficients are positive, which has to vanish. Thus we can conclude that 
the only possible solutions have 7 < 1, i.e. these are all unphysical. 

4.4- Vanishing magnetic part of the Weyl tensor 

Now we have to impose the condition a = k, so that u"^ + = (/x + 2p)/3. If we 
differentiate the last expression, we obtain yu(p) = p + /xq, where /xq is some integration 
constant. With p as our coordinate, the IC can then be integrated to give 

E{p) = {2p + /io)/3 , uo\p) = ujl{2p + /io) , a\p) = -/io/6 + D{2p + /io) . 

This model belongs to a class described by Cahen and DeFrise in Ref. |2l 

4-5. Vanishing k 

From (4.2) we have u'^ = (/i+p)/2, which, as in the preceding cases, can be differentiated, 
and it will yield ^{p) = —3p + /io, if a 7^ 0. Here /io is an integration constant. E and 
H are give by Eqs. (4.5) and (4.6). 

Choosing p as our coordinate, we can integrate Eq. (4.3) if we insert the equation 
of state and the expression for tu^. The result is 

a\p)=p + C/i-2p + fio), (4.11) 



where C is an integration constant. In Ref. [£2[ there remained a second order 
differential equation [Eq. (65)] to solve. This is eliminated by the above choice of 
coordinate. 



5. Spatial twist in ci: A 7^ 



Proceeding as in Case 1, we obtain the IC as the differential equations 

XP = X^-{fL + p)/2 + p{a-(3) , (5.1) 

Xa= - 2\^ + (/i - p)/2 - a{2p + a) , (5.2) 

XX = X{a - 2(3) , (5.3) 

X/i = - (2/? + a)(/i + p) , (5.4) 

where ' denotes differentiation with respect to and the algebraic relations 

E = - 2A2 + 2/i/3 - 2pa , (5.5) 

H = 2X{a-p). (5.6) 
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As in Sec. ^, we can reduce this system to three differential equations for three functions 
in the generic case. 

5.1. Constant A 

When A = 0, Eq. (5.3) implies a = 2(3. Inserting this into Eqs. (5.1) and (5.2) gives 

= -2AV5 + (3/i + p)/20 . (5.7) 

This equation can be differentiated to give [by using the IC and assuming p = p{fi)] 
either (i) /3 = or (ii) 2{fj, + p)p^^ — {f^ + 3p) + 12A^ = 0, where (i) gives a homogeneous 
model. 

As a simple ansatz for (ii) we might try p = A/i + B, which gives (a) p = /x + 12A^, 
or (b) p = -fi/2 + 3A2. 

5.2. a = c(3, c constant 

Within this class of models we have the shear-free and expansion-free cases. Inserting 
the ansatz into the IC gives us 

2(2c + l)c(3'^ + (2A2 -p- ii)c + A\^ +p~ ii = Q , (5.8) 
8(c + 1)A2 + (p,^ - l)(/i + p)c + 2(/i + p)p^^ - 4p] (c - 1) = . (5.9) 

We see that this equation is automatically satisfied if c = 1, i.e., if the fluid is shear 
free. Finally, inserting Eq. (5.9) into Eq. (5.3) gives us, with the help of the rest of the 
IC, the following second order differential equation for the equation of state: 

{(c + 2)2 (/i + p) [(^ + p)p^^^ + {p^^f 

+2(c - 3)(c + 2)(/i + p)p,^ - 3c2(/i + p) + 2c(/i - 3p) + 16^} 

X {(c + 2)2(;i + p)p^^ + 6c(;U-p)+4(/i-3p)}(c-l) =0 (5.10) 

Once again, the above equation is automatically satisfied if c = 1. Also, if c = —2 (i.e., 
the expansion free case), we get rid of all derivative terms in Eq. (5.10). The remaining 
equation to solve is (5.4). One may make the ansatz p = (7 — l)/i for the equation of 
state. From Eq. (5.10) we then express c in terms of 7. There are three possibilities: 
(i) 7(2 + 7)0^-1-2(272-57 + 6)0 + 4(2-7)2 = 0, (ii) c = 2(2-7)/(2 + 7), or (iii) c = 1. 
Below we investigate some special values of c. 

5.2.1. c = Since we have divided by c at several places, we here use Eqs. (5.1)-(5.4) 
directly. The algebraic expression for A, which we obtain from Eq. (5.2), is differentiated, 
which gives (/i + p)p^^ + /i — 3p = 0, where an equation of state is assumed. This has 
the solution 

2/i/ (/i — p) + In l/x — p| = constant . (5-11) 
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Choosing /i as the coordinate, we have 

where A is a constant and /(/x) = dfl/[fl + 

5.2.2. c—1, shear free case Equation (5.10) now becomes automatically zero. Prom 
Eqs. (5.5) and (5.6) we find that E — H — 0. Thus we have an H3 and a Gq, 
which implies that we are deahng with the open FLRW models. (Note that A is now 
redundant . ) 

Of the above the above equations, the only remaining are Eqs. (5.3) and (5.4), while 



(5.13) 



Since we have no constraint on the equation of state, we choose /i as our coordinate 
(assuming it is non-constant). This implies, through Eq. (5.4), that X — — 3/3(// + p). 
Equation (5.3) can now be integrated to yield 



(5.14) 



(the case P — is space-time homogeneous). Here A is some constant of integration. 
As an example, we can solve this equation for p = (7 — 1)^, which will give 



A 



Here we have used — 3(5. 



5.2.3. c = —2, expansion free case We have /i = 0, and it is then straightforward to 
show that this is a homogeneous model. The homogeneity of the model occurs because 
of our assumption of an equation of state. We might take another route by not assuming 
this. We then have 



Xp = 8(3 \{i^-p)/4- 



(5.15) 
(5.16) 



Equation (5.3) takes the form XX — — 4/3A. We choose A as our coordinate. Then we 
can integrate Eq. (5.16): 

p^li- 2X^/5 - CVX , (5.17) 

where C is some constant. All other quantities are then given in terms of A: 

E = -2X^ + 2/X/3 - (4// - 6AV5 - 3CVa) /3 , ^ 3 (4/x - 6AV5 - 3CVa) A^ . 
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5. 3. Vanishing electric part of the Weyl tensor 

We now have the constraint i? = 0-vv-A^ — /i/3 + (3a = 0. Differentiating E = and 
using Eq. (5.3), we get two cases: (i) P = a, i.e., no shear, or (ii) 5/i — p = ISPa. Since 
case (i) has already been treated above, we concentrate on case (ii). From now on we 
assume fi ^ 0. Using (ii) and the IC leads to the equations 

a = — (3 , (5.18) 

.2 ^ 1 (/i + 3p)(l-2g) 

^ 12 2(l + g)g,^ + Q(l-2Q) ' ^ ■ ^ 

where Q =^ (1 — p,;i)/(3 + p,^). But by using (5.18) in the constraint (ii) we get 
= (5/i — p)/{36Q). Thus we obtain a constraint equation for the equation of state: 

16(5/i - p)p,^^ + (3 + p,^) (p,^ - 1) (1 + 3p,,) (/i - 5p) = . (5.20) 

We can solve by making an ansatz on the equation of state. We choose p = (7 — 1)/^. 
Inserting this into Eq. (|5.20|) gives 7 = 6/5, and our equation of state becomes p = fi/5. 
From Eqs. (5.19) and (5.18) we get 

P"^ = 8/i/15 , = 2/i/15 , 

respectively. This solution was found by Collins and Stewart 



5.4- Vanishing magnetic part of the Weyl tensor 

From (5.6) we get (3 = a, which implies H = 0. This leads to X(3 = Xa =^ 
— yu/3 + = 0. But from (5.5) we see that E = 0, which implies that this is 
the shear free case. Thus these are the open FLRW models. 



6. Homogeneous space-times 

This class of models is determined by algebraic equations only. We find that a = 
2P + a = from Eq. (|3^ ). Inserting this into the algebraic equations we deduce that 
H = 0. Also, we obtain a number of very simple relations such as ujS = ujk = uX = 
Pn = PX = Xk = Ek = {3E + fi + p)P = 0, plus three equations determining the 
relation between the kinematic and geometric quantities. From these equations we get 
the solutions 



1. K = X = a = H = 0, and cu^ = 3-E'/2 = fi = p, which is Godels model |2^ . 



2. K = X = uj = H = 0, and 3/5^ = 3E/2 = —fi = —p, so that we have the restriction 
of stiff matter with negative pressure and energy density. 
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3. K = uj = P = H = E = 0, fi + 3p = 0, and = /i/3, which is Einstein's 
static model. Since this has a Gj, the rotation coefficient A is redundant in this 
description. 

4. X = uj = f3 = H = E = 0, fi + 3p = 0, and = p. This is once more Einstein's 
static modeL Here, the rotation coefficient k is redundant. 



6.1. Dependence on one null coordinate 

We can choose x° as our null coordinate, so that x = y = (since everything is symmetric 
w.r.t. and x^, it does not matter which one of them we choose). That x^ is null means 
that X = eY, where e = ±1. Next we insert this into the IC. This gives a = —ea and 
uj = eX. 

By using these equations in the IC it is straightforward to show that these space- 
times are homogeneous, uj = \ = (3 = a = E = H = fi = p = Q and /i = — 3j9 together 
with = p, i.e., this is Einstein's static model. 



7. The isometry algebra 

With the method described above it is not only possible to find the dynamics of the 
fluid, but it is also possible to find the metric in the specific cases of interest. One uses 
the structure of the isometry group on the orbits in F{M) to make a metric ansatz. 
This is an easy task in the case of three dimensional groups, since one then uses the 
Bianchi classification, and it also works well in the four dimensional case. From this one 
can calculate the reduced set S and compare with the original set and from this obtain 
the metric components in terms of the kinematic quantities, fluid variables and rotation 
coefficients. 

We start out with Cartan's equations (2.3) and (2.4), and project these onto 
the orbits of the isometry group in F{M). Working in F{M) has the advantage of 
always giving an isometry group acting simply transitive This in turn means that 
the differential algebra between the basis 1-forms (Cartan's equations) generates the 
structure constants of the isometry group (if we choose an invariant basis {e/}), i.e., if 
{^p} are the Killing vectors in F{M) then 

Ip,Iq] = C%^n ^ du^'l = ^-C'^RsOJ^'l A , (7.1) 

where u}^\ is the projection onto the cotangent space of the orbits in F[M) and the 
C^RS are defined in terms of C^jk through 

C^RS = C^Rs + C^abI\I^\rI\I^\s — 2C^Asl\l^\R ■ (7.2) 

Thus the rotation coefficients and the Riemann tensor [on the orbits in F{M)] essentially 
corresponds to the structure constants of the isometry group of the space-time. 
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The orbits in F{M) are defined as d/^ = 0, or when using a fixed frame dx" = 0, 
T% = 0. This means that some of the 1-forms will be linearly dependent on the others. 
In the case of LRS, and with our choice of frame, this means that, in the generic case, 

is the connection 1-form that it linearly independent of the cotangent basis. 

In the following sections we find the metrics for the first two LRS classes, and discuss 
their relations to the IC. 

8. The isometry algebra for cases 1 and 2 

The isometry algebra [Eq. (7.1)] (when 7^ or A 7^ 0) can be written as 



d(T^ = T(T^ A(T^ , (8.1) 

do-2 = a^A(T\ (8.2) 

do-3 = a^Aa\ (8.3) 

do-^ = Scr2 A , (8.4) 



where cr = u \, A = 2,3, and the definitions of the remaining quantities can be 
found in Table 1. 

Table 1. The definitions of the quantities used in Sees. ||and||. Notice that the main 
diff'erence between the two cases hes in the causal properties of the essential coordinate 



Case 






cr 




r 


S 








q 


V 


1. CJ 7^ 






















time-like 


2. A 7^ 








a;" 


2A 






+ 3A2 






space- like 



We may regard Eq. (8.1) as a separate (differential) equation for cr^, and Eqs. 
(8.2)-(8.4) as the algebra we need to solve. can then be found through cr^. When 
7^ 0, we know that ux" = y^^dx^, and when A 7^ we know that uj^ = X~^dx^. 
The isometry algebra can be split into three cases, according to whether E > 0, = 
or < 0. 



8.1. S^O 

If we 'normalize' our 1-forms as = -\/Scr^, Eqs. (8.2)-(8.4) obtain the structure of 

Bianchi type IX. The canonical solution is (see Ref. |TT| ) 

cr^ = — sin (^d6' -|- sin 6* cos Cd0 , (8.5) 

cr^ = cos Cd6' sin 61 sin Cd0 , (8.6) 

(7^ = cosM0 + dC , (8.7) 
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where 9, 0, and ( are some coordinates. The equation for cr^ becomes dcr^ — 
(r/E)cr^ A cr^. With the result above, we can solve for <t^ and obtain 

<T^ = (r/S)cosM(/. , (8.8) 

where we have neglected a total differential in the integration procedure. We may 
proceed in two ways according to the properties of E. 

If E > 0, the coordinates in the solution are real valued, and we can apply it directly. 

If E < 0, the new 1-forms introduced above take their values over the set of complex 
functions. Thus we need to make the transformation 9 — > i9, where the 'new' 9 is real. 
This will introduce an overall i in the solutions for the 1-forms (since sin ^ — > i sinh 9) , 
which cancels the i in the definition of the 'normalized' 1-forms. Also, in the solution 
for <T^, cos^ cosh 6', due to the transformation. 



8.2. E = 

We can now define <t^ <T^±i<T'^ and cr^ —icr^. Equations (8.2)-(8.4) then becomes 

d<T+ = (T+ A (T^ , (8.9) 

dcr~ = -cr~ Act'' , (8.10) 

d<T^ = . (8.11) 

This has the structure of Bianchi type VI, and has the solution 

0-+ = e-^dri , = e^dr , o"^ = d^ . (8.12) 

Since the original 1-forms takes their values over the real functions, we make the 

transformation ^ = — i^, r] — t* = 9 + i(f), for some real coordinates 9, and (p 
(here * denotes complex conjugation). Thus we can write 

= cos Cd^ - sin Cd0 , (8.13) 

= sin Cd^ + cos Cd0 , (8.14) 

o-^ = dC . (8.15) 

The equation for cr^ becomes da^ = {i/2)Tcr'^ A <t~, which has the solution 

(T^ = r9d(l) . (8.16) 

9. Metrics for cases 1 and 2 

We here make the ansatze for the full metric for space-times where a; or A 7^ 0. 
When T, ^ the simplest guess is to take 

(T^j:x{qW + tp{q)dv , (9.1) 

= S{q)&^ . (9.2) 
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When E = 0, we can make an ansatz analogous to the one above: 

<j = x{q)<T'+il^{.q)dv , (9.3) 
= 5{q)cr^ . (9.4) 

(See Table 1 for the definitions of quantities introduced above.) 
9.1. Metrics 

Introducing the parameter e, which is defined as —1 for case 1, and 1 for case 2, we can 
write the metric g = e(crc)^ ~*"^''g) where <Tc is the 'complement' to <j (see Table 1), in 
one common form for all the cases. The metric on surfaces {dg = 0} becomes 



(3)o 



' e {il^dv + xr cos ed(t)f + 52 (d^2 ^ gjn2 ^^^2^ ^ Y.>Q 

e (^dt; + xrM0)^ + 52 (d^2 _^ ^02) ^ ^ q 

_ e (V^dz; + xr cosh ^d0)^ + 5"^ (d^^ gj^^^ ^dc/)^) , E < 



d^^ + 4-^' (e^^-sV^^)'d(/.^ 



(9.5) 



Here we have introduced another parameter, s, which is defined as —1 when E > 0, 
when E = 0, and 1 when E < 0. 

9.2. Consistency equations 

With consistency equations we mean a comparison between the old set S and the new 
set 5* calculated from the metric ansatz. As an explicit example of how to do this, we 
will look at case 1. 

Starting from our metric ansatz, we use the type of tetrad that lead to the ansatz (in 
this case a Lorentz-tetrad) and calculate the rotation coefficients T%i. We then observe 
that ro23 = a;x/5^, which then gives x = by comparison. Inserting this into Foia 
and comparing with Foia from 5", we get the relation cuS^/ip — constant. This constant 
can be absorbed by redefining our time coordinate by a simple scaling (and refiection). 
Doing this we obtain i/j = uj5^, i.e., the only function from the ansatz that remains to 
be determined is S (the others are known from the IC or their relation to 5). 

For case 1, the remaining equations are 



n = -y5' /5, 

E-{pi + 3p)/6 = -y [{yco'y/u; + 2{y5y/5] - Ay'5'u;' /{5u;) - 2{y5'/5)\ 



(9.6) 
(9.7) 
(9.8) 
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-E/2 - (// + 3p) /6 = -u;' - y^S'u'/iSu;) - 2{y5'/5f, (9.9) 
£;/2-/i/3 = y(y5')'A (9.10) 
- ///3 = -So;' + + (|/575)2, (9.11) 
= -2t/(5cj)V5 , (9.12) 
where s was introduced in the general metric (9.5). 



If one insert the expressions for a and k into the remaining consistency equations, 
and then use the IC, one discovers that the only equation that is not trivially satisfied 
when s 7^ is Eq. (9.11). Thus this equation determines the metric when s 7^ 0, ones 
the IC are solved. If s = 0, Eq. (9.7) determines 5. Case 2 can be treated in the same 
manner, and we give the function 5 for the different cases in Table 2. 



Table 2. The form of the function 5 for the different cases. Here A is an integration 

constant. 



Case 


S 7^ 




E 


= 


I. u^O 


52 = s{-E- 


Ai/3H 


-3co2-k2)-i = >0 S 


= A cxp [— J (K/j/)dx^] 


2. A ^ 


6^ = s{-E- 


/x/3H 


-3A2 + /32)-i = |5]|-i > 6 


= Aexp [j\(3/X)Ax°] 



The metric (9.5) can now be written as 



+ 5^ 



d^' + 4- 



(9.13) 



where the function 5 is given in Table 2. Having solved the IC, we know the metric 
explicitly, since all functions are determined, including 5. 



10. The function S and its relation to the IC 



Since we have three different metrics depending on the properties of the ftmction E, 
we must, through the IC, get some restriction on the constants of integration that has 
been introduced when solving the IC. When sgn(S) = ±1 we will get some inequality 
that have to be fulfilled by the integration constants. This has to be treated separately 
for every solution of the IC. On the other hand, when S = 0, we have an algebraic 
constraint on the kinematic quantities. This will lead to more severe restrictions than 
the inequalities impose. We can therefore treat this case on its own. 

10.1. Case l,ujy^O 

Now we can write — T, — E + iJ,/3 — + — —p + 2aK — + k^. Taking the 
derivative of this and using the IC, gives us either (i) a = 0, or (ii) k — 0. 
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Case (i) was treated in Sec. but now there is the extra equation uj^ = —p + h?. 
The energy density must not be constant, since then = 0, i.e. a contradiction to our 
starting assumption. If the energy density is non-constant, we choose it as coordinate. 
We then obtain the solution of Sec. ^.2| , but because of our extra constraint S = 0, we 
find that the integration constant A = 0. 

For (ii), we obtain u"^ = —p, so that the pressure is negative. We find that p must 
not be constant, because then u = 0. Thus this case is equivalent to the one treated in 
Sec. with the constant of integration /iq = 0. 



10.2. Case 2, X ^ 

The constraint becomes = S = E+^/S + SX^ — P"^ = A^+/i — /9(/3+2a). Differentiating 
this and using the IC gives an identity. Thus we can get rid of, say, A (expressing it 
algebraically in terms of the other functions) and get a simpler set of integrability 
conditions. The IC may then be solved (as before) for some special cases, e.g. an 
expansion free fluid. 



11. Concluding remarks 

As a conclusion one might point on the advantages of this method compared with other 
methods. First, we automatically obtain a reduced set of equations (although they 
may contain redundant information) needed to be solved to find a solution to Einstein's 
equations. Second, we have a split in the equations, one set consisting of the equations 
mentioned above, and one set that determines the rest of the metric [the 'redundant' 
part connected to the symmetry of (M, g)]. We hope that these examples have shown 
on some of the strengths of the method presented. 
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Appendix A. Integrability conditions 



In our fixed frame formalism we get the following IC: 
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1. Commutator equations x°'\[k^\i3\x^\ij — x°'\mr"^[ki], with the indices (^J given in the 
bracket: 



XY + X'y -YX- Y'x ^Xa^Ya , 
(23) = Xa;-rA, 
(^01) xY + x!y — yX — y'x = + ya, 
(23) Q^xu- yX, 



(Al) 
(A2) 
(A3) 
(A4) 



2. Ricci identities = 2r'^6[j- + 2r'''"[jr|b„|i] + 2r"fefcr^'[ij], with the indices 

(^jj) given in the bracket: 



\ 

101 J 


Ya + ya' -- 


= + - (// + 3p)/6 - - - xa' , 


(A5) 


\ 
123 J 


H^2(a- 


- k)uj + 2{a - (3)X, 


(A6) 


\ 
202 J 


Xl3 + xl3' 


= -{fi + 3p)/6- - E/2 + aK + u;^, 


(A7) 


\ 
212 J 


YP + yP' 


= k{P — a) + Xu, 


(A8) 


\ 
203 J 


Xuj + xuj' 


= -aX - 2u;p, 


(A9) 


A 

213 J 


Yuj + yuj' 


= -H/2 + KW + X{a-P), 


(AlO) 


202) 


Xk + xk' 


— {a — k)P — Xcu, 


(All) 


212) 


Yk + yn' 


= -E/2 + - A^ + /i/3 - Pa, 


(A12) 


203) 


XX + xA' 


= -{a- K)uj + H/2-Xp, 


(A13) 


213) 


YX + yA' 


— 2kX + cua, 


(A14) 



3. Bianchi identities R^q[ij;k] 
bracket: 



0, with {''g} = {%}, and the indices djk) given in the 



(023) x(^E + fi/3)+x{E' + i2'/3) = -{3E + ii + p)p -3HX, 
(123) y[e + ii/i) + y{E' + At73) = ?>Ek - SHcu. 



(A15) 
(A16) 



Here we have used the notation ' = d/dx^ and ' = d/dx^ . 

We observe that the cychc identities R^[ijk] = are aheady imposed by our choice 
of Riemann tensor. 
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